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Abstract 

A new type of exact solvability is reported. Schrodinger equation is considered in a 
very large spatial dimension D ^> 1 and its central polynomial potential is allowed 
to depend on "many" (= 2q) coupling constants. In a search for its bound states 
possessing an exact and elementary wave functions ip (proportional to a harmonic- 
oscillator-like polynomial of a freely varying, i.e., not just small, degree N), the 
"solvability conditions" are known to form a complicated nonlinear set which requires 
a purely numerical treatment at a generic choice of D, q and N. Assuming that D is 
large we discovered and demonstrate that this problem may be completely factorized 
and acquires an amazingly simple exact solution at all N and up to q = 5 at least. 

PACS 03.65. Ge 



1 Introduction 



The key motivation of our present study of polynomial oscillators lies in the well 
known fact that the majority of quantitative predictions in nuclear, atomic, molec- 
ular and condensed matter physics must rely on a more or less purely numerical 
model. The completely non-numerically tractable quantum systems are rare though, 
at the same time, useful and transparent (cf., e.g., the description of vibrations in 
molecules mimicked by harmonic oscillator). Polynomial oscillators may be taken, 
in this setting, as lying somewhere on a borderline between the two regimes. 

The first indications of a breakdown of the traditional separation between the 
numerical and analytic models in quantum mechanics came with the emergence of 
certain "incompletely analytic" (now we call them quasi-exactly solvable, QES) poly- 
nomial oscillator models. Unfortunately, their separate discoveries at the beginning 
of the last quarter of the twentieth century [1] - [5] have all been understood and 
accepted as a mere formal curiosity. As an interesting confirmation of analogies with 
classical mechanics for the charged and shifted oscillator in two and three dimen- 
sions [1], as a peculiar singularity in a universal continued-fraction algorithm for 
sextic oscillators [2] , as an exceptional case in the standard infinite-series solution of 
the corresponding differential Schrodinger equations [3] - [5] etc. 

The defining property of the QES models (being solvable just for a part of the 
complete set of their coupling constants and/or energies) re-acquired a new mean- 
ing only after their Lie-algebraic re-interpretation [6] which revealed their genuine 
mathematical closeness to the exactly solvable (ES) Hamiltonians [7]. Subsequently, 
in physics, a multitude of their relationship to the apparently different models has 
been revealed in a way summarized, e.g., in the detailed Ushveridze's monograph [8]. 

In spite of the unique success of the mathematical QES formulae, a number of 
difficulties remained connected with their practical applications and applicability. 
One of the key reasons (and differences from the current ES models) is that the 
explicit construction of the multiplet QES energies remains purely numerical. Indeed, 
these values must be computed as roots of an N— dimensional secular determinant so 
that the difference between the variational, "generic N = oo" rule in Hilbert space 
seems only marginally simplified by the QES construction of any QES multiplet with 

The main purpose of our present study is related precisely to the latter point. Our 
idea may be explained, briefly, as an application of perturbative philosophy assuming 
that the spatial dimension D is large. In this spirit we are going to consider a general 
Magyari's [4] QES Hamiltonian H^ N \D) (at a fixed dimension D of the space of 
coordinates and with a chosen size N of its secular determinants, see below for a 
more detailed explanation). Finally we construct a set of its specific approximations 
H^ N \oo) with errors proportional to 1/D. 

The exact solvability of the latter Hamiltonians Hq 1 '^ emerged as an utterly 
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unexpected result of our (originally, fully numerical) calculations. Our presentation 
of it starts from a concise review of the concepts of exact solvability in section 2. 
Section 3 will then mention a few specific formal merits of transition to the domain of 
the large dimensions D which simplifies the general Magyari's secular determinants 
considerably. The core of our message appears in section 4 where in the domain of 
D >> 1, the separate families of the polynomial oscillators (numbered by the integer 
q — 1,2 and 3) are studied in detail and shown to lead to the closed solutions. 
Section 5 outlines the possibilities of an extension of these results to q = 4 and 5 
while section 6 summarizes and discusses some possible practical consequences. 



2 Exactly solvable oscillators: A brief review 

2.1 Harmonic oscillator and the like: All bound states are 
elementary 

One of the most exceptional exactly solvable models in quantum mechanics is the 
central harmonic oscillator in D dimensions. Its so called superintegrability (the term 
coined by Pavel Winternitz [9]) makes its partial differential Schrodinger equation 

(-^-/\ + ]-mn 2 \x\ 2 \m(x) (1) 

solvable by the separation of variables in several systems of coordinates. The most 
common cartesian choice may be recommended for the first few lowest spatial dimen- 
sions D only [10]. In contrast, the separation in spherical system remains equally 
transparent at any D because it reduces eq. (1) to the same ordinary (so called 
radial) differential equation 



d 2 1(1 +1) 2 2 

h — + urr 

dr 2 r 2 



ip(r) — Eip(r) (2) 



with r = \x\ E (0, oo), E = 2me/h 2 and oo = mVL/h > 0. In this language we have 
I = £ L = L + (D — 3)/2 where L — 0, 1, . . . [11]. At each L the energy levels are 
numbered by the second integer, 

E = E n>L =u(2n + e L + 3/2), n, L = 0, 1, . . . . (3) 

The wave functions with quadratic A(r) = ur 2 /2 > and minimal iV = n + 1 in 

AT-l 

Mr) = r e+1 e- A « £ h m r 2m (4) 

m=0 
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are proportional to an n— th Laguerre polynomial [12]. In Hilbert space, their set 
is complete: This may be explained via oscillation theorems and characterizes the 
harmonic oscillator as exceptional. In such a setting, one should recollect all the sim- 
ilar (i.e., Coulomb and Morse) exactly solvable potentials, but one need not mention 
them separately as long as they are formally equivalent to our harmonic oscillator 
example after a simple change of variables [13]. 

2.2 Sextic oscillators: Many bound states are elementary 

For many phenomenological purposes, even the shifted harmonic oscillator forces 

+ uj 2 r 2 are not flexible enough. Fortunately, there exists an 
immediate QES generalization of the harmonic oscillators revealed by Singh et al [2] . 
Let us summarize briefly this type of solvability as a construction which starts from 
a replacement of the free constants in V^ HO \r) by their polynomial descendants of 
the first order in r 2 , 

uj — > W(r) = a + a x r 2 , V {HO) (0) — > G_i + G r 2 = U(r) . (5) 

After the conventional choice of G-i = this trick gives the general sextic potential 

V (sextic \r) = U(r) +r 2 W 2 (r) = g r 2 + 9l r 4 + g 2 r 6 (6) 

all three couplings of which are simple functions of our initial three parameters and 
vice versa, 

92 = ct\ > 0, g 1 = 2a a 1 , g = 2al + G . (7) 

The resulting Schrodinger bound state problem cannot be solved in closed form. 
Nevertheless, we may postulate the polynomiality of the wave functions ip^T^ ( r ) f or 
a finite multiplet (i.e., iV— plet) of the wave functions. For this purpose it is necessary 
to lower the number of the freely variable couplings by the specific constraint 

G = -a 2 - ai (4N + 2C+1), N>1. (8) 

For eachiV— plet, the polynomial solutions (4) are made exact by the N— dependent 
QES condition (8). The choice of a WKB-like (i.e., quartic) exponent 

1 1 

A(r) = ^"o^ 2 + 4«i r4 ( 9 ) 

guarantees their physical normalizability. The ansatz (4) transforms then the differ- 
ential Schrodinger equation into a linear algebraic definition of the unknown iV— plet 
of coefficients h m . An incomplete solution is always obtained for a mere finite set 
of the levels n G (n ,ni, . . . , njv-i)- In contrast to the harmonic oscillator, the QES 
solvability is based on the L— and N— dependent constraint (8) so that, generically, 
the elementary QES multiplet exists in a single partial wave only. 
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2.3 General harmonic-oscillator-like bound states 

The explicit energy formula (3) for harmonic oscillator was replaced by an implicit 
definition in the preceding paragraph which gives the sextic QES energies in the 
purely numerical form, viz., as zeros of the Singh's secular determinant of a certain 
tridiagonal N by N matrix [8]. In this sense, a very natural further extension of 
the Singh's QES construction exists and has been described by Magyari [4]. Its 
description may start from the more consequent replacement (5), constructed from 
even polynomials of degree 2q, 

V^\r) = U {q \r) + r 2 [W {q \r)} 2 , U {q \r) = G r 2 + dr 4 + . . . + G q ^r 2q , 

W [q \r) = a G + a l r 2 + ... + a q r 2q (10) 
This formula re-parametrizes the polynomial 

V® (r)=g r 2 + g 1 r i + ... + g 2q r iq+2 g 2q = 7 2 > (11) 

and specifies the one-to-one correspondence between the two sets of couplings, 

{go, ■ ■ ■ , 92q} "<=^ {Go, • • • , Gq-1, «0, • • • , Ciq] 

where g 2q = a q 2 , g 2q -i = g2q-i{ot q , a. q -i) = 2« r i a„... or, in opposite direction, 
oiq = ^[gTq = 7 > 0, a,_i = g 2q -i/ (2a,) etc. 

At a generic q = 1,2,..., equation (9) must be further modified in such a way 
that rW{r) = Z'(r), 

A«(r) = \a Q r 2 + + ... + ^^a q r 2q+2 . (12) 

With a q > 0, one verifies that 

^(j>hysical)^^ e -XM(r) + 0(l)^ r>1 

which means that the correct bound-state ansatz 

N-l 

^(r) = J2 hl N) r 2n+e+1 exp {-\ {q) (r)] (13) 

n=0 

converts our radial equation (2) + (11) into an equivalent linear algebraic problem 

Q^h^ = 0. (14) 
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Its closer inspection reveals that this problem is overcomplete, i.e., its matrix is 
non-square and asymmetric, 



( B C 
A? B x d 



Q [N] 



q 



1(9) 

l N-2 



B n 



C n 



A 



(9) 

N-l 



4W R 
S± N _2 -DJV-2 



A 



(i) 

N-l 



Cn-2 
Bn-i 



A 



(9) 

N+q-2 ^N+q-2 



A 



(9-1) 



V 



A 



(9) 

N+q-l 



(15) 



Its elements depend on the parameters in bilinear manner, 



C n = (2n + 2) (2n + 2£ + 3), B n = E - a (An + 21 + 3) 
A« = -«i (4n + 2£ + 1) + al - g , A® = -a 2 (4n + 2£ - 1) + 2a «i - ^1, 

• ' ' 1 

A$ = —a q (An + 21 + 3 - 2q) + (a Q; 9 _i + aia q - 2 + • • • + a q -ioto) - g q -i, 

71 = 0,1,.... 

(16) 

At any fixed and finite dimension N = 1,2,... the non-square system (14) is an 
over-determined set of N + q linear equations for the iV non-vanishing components 
of the vector h^ N \ 



2.4 Changes of variables and an extension of applicability of 
the harmonic-oscillator-like constructions 

At q — 0, equations (15) degenerate back to the recurrences and define the well 
known harmonic oscillator states. As already mentioned, their additional merit lies 
in the coincidence of their polynomial part with the current Laguere polynomials. 

At q = 1 we return to the sextic model where the "redundant" last row fixes 
one of the couplings and where we are left with a diagonalization of an iV by N 
matrix which defines the N— plet of the real QES energies in principle. In such a 
setting, an important piece of an additional encouragement results from the well 
known possibility of a definition of new QES Hamiltonians by a change of variables 
r — > x and ip(r) — > x const 4>(x) in the Schrodinger equation [13, 14]. Even when 
considered just in the most elementary power-law form, this change is defined by the 
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prescription [15] 



r 2j _^ x 5(k)^ 6 ^ = 2 l±l- 2 , fc = l,2,...,2g + 2, 

K 

and extends the class of the " interesting" potentials not only in the well known q = 
case (see the open possibility of a transition to the completely solvable Coulombic 
potential as mentioned above) but also in the q — 1 model where the following 
four equivalent potentials may be distinguished and numbered by the above index k 
attached to them as their second superscript, 

y(5=l,fc=l)( r ) = ar 2 + &r 4 + r 6 
y(q=l,k=2)^ j = aj .-l + br + r 2 
y(q=l,k=3) ( r ) = ar -4/3 + fo.-2/3 + r 2/3 U 7 ) 

y(«=i.*=4)( r ) = ar~ 3 / 2 + fer" 1 + cr" 1 / 2 . 

The situation is more complicated at g > 1. The counting of parameters and equa- 
tions indicates that unless one broadens the class of potentials, only a very small 
multiplet of bound states may remain available in closed form [14]. Still, the same 
elementary change of variables enables us to extend the set of the partially solvable 
potentials to the six- member family at q — 2, 



V( q=2 > k=1 \r) = ar 2 + br 4 + cr 6 + dr 8 + 
V iq=2 > k=2 \r) = ar' 1 + br + cr 2 + dr 3 + . 

y( q =2,k=3) (r) = ar -4/3 + br -2/3 + cr 2/3 + rfr 4/3 + f .2 
y( q =2,k=4)( r } = ar -3/2 + 6r ~l + cr -l/2 + ^1/2 + r l 1 ^ 
y(q=2,k=5)^ = ar S/5 + &r -6/5 + cr -4/5 + ^-2/5 + r 2/5 
y( 9 =2,fe=6) (r) = ar -5/3 + 6r -4/3 + cr -l + rfr -2/3 + f r ~l/3 

etc (cf., e.g., paper [15] where the similar lists of the partially solvable potentials have 
been displayed up to q = 4). In this way, the availability of the exact solutions for 
all these various forces might offer a new inspiration, say, in some phenomenological 
considerations and models and/or for their perturbative analyses and some more 
detailed large— D calculations. 



3 Solvability of polynomial oscillators at the large 
spatial dimensions 

Up to now, our attention has been concentrated upon the structure of the QES wave 
functions. From the point of view of the evaluation of the energies, the main dividing 
line between the solvable and unsolvable spectra is in fact marked by the distinction 
between the closed q = formulae and their implicit QES form at q = 1. The 
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transition to the next q — 2 may be perceived as merely technical. At all q > 1, 
the difficulties grow with N. In such a setting we noticed the emergence of certain 
simplifications at D > 1. 

3.1 Difficulties arising at q > 1 

At any D, the last row in eq. (14) decouples from the rest of the system. At any 
q > 1 it may treated as a constraint which generalizes eq. (8), 

g q -i = —a q (An + 21 + 3 - 2q) + (a a q _i + a x a q ^ 2 + • • • + a^ao) . (19) 

The insertion of this explicit definition of the coupling g q _i simplifies the lowest 
diagonal in Q^ N \ 

A® =47(JV + g-ra-l). (20) 

Since A$ + x = we may omit the last line from eq. (15) and drop the "hat" " of 
the diminished matrix Q^ N K This gives eq. (14) in the more compact form 

QW fr N) = (21) 

where the size of the non-square matrix is merely (N + q — 1) by N. Unfortu- 
nately, the new equation is still purely numerical, with an exception of the simplest 
special case with q = where no coupling is fixed and where the energies themselves 
are given by the explicit formula (19). At q = also the recurrences for coefficients 
of the wave functions may be solved in a compact form. 

As already mentioned, the q — 1 version of eq. (21) degenerates to the secular 
equation 

det Q [N] = 0. (22) 

This is a purely numerical problem at all the larger N > 5. Still, one coupling is 
fixed by eq. (19) and only the iV— plet of energies must be calculated as the real 
zeros of secular polynomial. 

At q > 2 the q independent (and mutually coupled) N by N secular determinants 
must vanish at once [16]. With an auxiliary abbreviation for the energy E = —g~\ 
this means that at least one of the couplings is always energy- dependent and its 
value must be determined numerically. In the other words, our non-square matrix 
QW = Q [N] (g -ij <7oj • • • ,9q-2) w iH annihilate the vector if and only if all its q 
arguments are determined in a deeply nonlinear and self-consistent, mostly purely 
numerical manner. 

3.2 Simplifications arising at D > 1 

For a clear understanding of what happens at D ^> 1, let us pick up the q = model 
V^(r) and re-consider its coordinate-dependence in the D > 1 regime. We discover 
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a quick growth of the minimum of the effective potential, occuring at a fairly large 
value of the coordinate R = R{D) = [1(1 + l)/u 2 } l/A » 1. Its Taylor expansion 

^-ti + u 2 r 2 = 2u 2 R 2 + 4cu 2 (r - R) 2 - l( r - i?) 3 + . . . (23) 

reveals that the shape of the effective potential is R— and D — independent. Near 
the minimum, also the cubic and higher corrections become negligible. This implies 
that the shifted harmonic oscillations characterize now the local solutions very well. 
In particular we have the wave functions 

^ ~ e — (r-V 2 ; fa~( r -R) e~" {r ~ R)2 , . . . , r^R (24) 

and re-derive also the leading-order degeneracy of the spectrum and its equidistance 
in the next order, 

E Q = 2u 2 R 2 + 2u + ..., E x = 2u 2 R 2 + 6u + . . . , .... (25) 

The agreement of this approximate formula with the available exact spectrum (3) is 
amazing. 

Considerations which have led to this agreement may be applied as a guide to the 
large-.D description of the low-lying states in any phenomenological input potential 
well V(r). In such a context, a skeptical question is due. Once we have the exact 
formula, why should we search for its alternative (re-) derivation? The reply will 
follow from our forthcoming results. We shall see that a qualitatively different new 
source of D ^> 1 simplification will emerge in all the exactly solvable polynomial 
q < oo models. 



3.3 Magyari equations in the large— D regime 

In the above D ^> 1 construction, little information can be extracted from the wave 
functions themselves. Although we Taylor-expanded the effective potential, we did 
not make any use of the information about the wave functions. In this sense, we 
are now going to demonstrate the feasibility of the approach where the guaranteed 
polynomiality of the wave functions will play a key role. 

In our original differential eq. (2) as well as in all its q > generalizations, 
the numerical value of the spatial dimension D will be assumed large. This will 
simplify our matrix re- arrangement (21) of this problem with the matrix elements 
(16) re- written as the linear functions of D, 

C n = (2n + 2) (2n + 2L + D), B n = -g^ - a (4n + 2L + D), 
A { n ] = -g k -i - a k (4n + 2L + D - 2k) + (a «fc-i + • • • + a fc _ia ) , (26) 
k = 1,2,. 1, n = 0,l,...,N + q-2 . 
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Besides the D— independent Affl = Affl^ which remains unchanged, we shall pre- 
serve the dominant components of the matrix elements, 

Cf = (2n + 2) D, B® = -<?_i - a D, = -g k ^ -a k D,k<q. 

Then we re-scale the coordinates and, hence, coefficients, 

= Pn /V n ■ (27) 

By the choice of the parameter /i we are free to achieve that the uppermost and the 
lowest diagonals are just a re-ordering of each other, 



2D 



r = 4 7 /A 



(28) 



This exhausts the freedom and fixes the D — dependence of our scaling, 



H = ti{D) 



,27) 



1/(9+1) 



t = t(D) = (2 9+2 £> 9 7 ) 



1/(9+1) 



(29) 



The recipe replaces the energies and couplings {g~i, go, ... , g q -z} by the new param- 
eters {si, S2, ■ ■ ■ , s q } in linear way, 



9k-2 = -afe-i-D - 



fc-i 



Sk, 



k= 1,2,. 



(30) 



In the leading-order approximation, our re-scaled Magyari equations read 



«2 



N - 1 



1 



\ 



N-2 s n 



V 



si N-2 

si N-l 

: si 

2 Sq Sq—1 

1 S n 



( po ^ 

pi 

PN-2 

\ Pn-i J 



= 



(31) 



We shall now study their solutions. 



3.4 Inspiration: New closed formulae at q = 

At q — 0, as already mentioned, the energies are unique functions of iV and do not 
exhibit any unexpected behaviour. Still, it is instructive to extract the leading-order 
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result for wave functions. The q = version of eq. (31) 



/ N- 1 



1 

N -2 2 



V 



N -2 

1 TV- 1 



Pi 

PN-2 



defines the (up to the normalization, unique) Taylor coefficients 

' N — 1 \ 



Pn=(-l) n pO 



n 



(32) 



(33) 



We may appreciate that they do not carry any round-off error and that the related 
leading-order wave functions possess the elementary form 



1 - 



-Xr 2 



n 



0,1, 



(34) 



At the same time one must be aware that only the leading-order part of (34) is to 
be compared with the available exact D < oo result. In particular, the presence 
of a degenerate nodal zero is an artifact of the zero-order construction. All this 
experience may serve as a guide to the less transparent q > cases. 



4 New partially solvable models with q > 1, any 
N and large D ^> 1 

In the way inspired by eq. (18), one may move beyond q = and q = 1 and transform 
the decadic forces into their quartic equivalents etc. Paper [15] may be consulted for 
details which indicate that the study of any potential V(r) which is a polynomial in 
the powers of the coordinate r may be replaced by the study of its present Magyari's 
or "canonical" QES representation V^ q \r) at a suitable integer q. In addition, we 
shall also restrict our attention to the domain of large D. 



4.1 Guide: Sextic QES oscillator with q = 1 and any N 

Starting from the first nontrivial sextic-oscillator potential (6) with q = 1 and with 
the binding energies re-parametrized in accord with eq. (30) where S\ = s, 

E= 1 -^D + (64g 2 ) 1 /*VDs, 
1 V92 
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full attention must be paid to the selfconsistency problem represented by the set of 
equations (31). At every N, its first nontrivial q — 1 version 



/ s 1 
N-l s 

N-2 



V 



\ 



2 

s 



3 
2 



N-l 

s 



( 



Po 

Pi 



\ 



PN-2 

V Pn-i J 



(35) 



has the form of an asymmetric eigenvalue problem. In standard manner it leads 
to the secular equation (22) expressible as the following sequence of the polynomial 
conditions, 

s 3 -4s = 0, N = 3, 

s 4 - 10 s 2 + 9 = 0, N = A, 

s 5 -20s 3 + 64s = 0, N = 5, 

etc. By mathematical induction, all the infinite hierarchy of these equations has been 
recently derived and solved in ref. [17]. 

Quite remarkably, all of the real (i.e., "physical") energy roots s = proved 
to be equal to integers. Moreover, all of them may be determined by the single and 
compact formula 



,(?) = _iv - 1 + 2j, 



J = 1,2, 



N. 



(36) 



This represented one of the key motivations of our present work, especially when we 
imagined that also all the related coefficients p$ may equally easily be normalized 
to integers, 

1, N = l, 



P ( l ] 



Pt ] = -P? = 1, 



N = 2, 



(1) (1) (2) (2) (3) 

PO = P2 = PO = -P2 = PO 

etc. 



V? = 1, 



(3) 
-Pi 



2, pS 2) =0, iV = 3, 



The first result of our subsequent computations using the symbolic manipulation 
techniques proved equally encouraging since we succeeded in compactification of the 
set of the above recurrent solutions to the single leading-order form of the related 
wave functions, 



i + 



r 
/i 



exp 



1 



1 



-a r 



-air 



j = l,2,. ..,N. 



(37) 



A few more comments may be added. 
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• The large and degenerate nodal zeros in eq. (37) are a mere artifact of the 
zero-order construction. 

- The apparently interesting exact summability of all the separate C(r 2 //i) 
error terms is not too relevant, indeed. Although it leads to the zero- 
order nodes at r = O(yjp) = 0{D 1 ^)^ these nodes have no real physical 
meaning. 

• The leading-order perturbative approximation provides a reliable information 
about the energies. 

- They are asymptotically degenerate, due to the large overall shift of the 
energy scale as explained in section 3.2. 

— The next-order corrections may be easily obtained by the recipes of the 
textbook perturbation theory. 

— As long as the coefficients p n are defined in integer arithmetics, the latter 
strategy gives, by construction, all the above-mentioned energy corrections 
without any rounding errors in a way outlined in more detail in ref. [17]. 

In the other words, we may say that formula (37) may either be truncated to its 
leading-order form jp^(r) = r e+1 exp (-A^(r)) or, better, its full form may be 
used as a generating function which facilitates the explicit evaluation of the coef- 
ficients p(i\ In comparison with the oversimplified harmonic oscillator, the q — 1 
wave functions may be characterized by the similar coordinate dependence which 
becomes spurious (i.e., dependent on the selected normalization) everywhere beyond 
the perturbatively accessible domain of r. 

At the same time, the energies specified by eq. (36) form and absolutely amaz- 
ing multiplet. On the background of its existence, a natural question arises whether 
some similar regularities could also emerge at some of the larger integer indices q > 1. 
We are now going to demonstrate that in spite of the growth of the technical obsta- 
cles in dealing with the corresponding key equation (31), the answer is, definitely, 
affirmative. 

4.2 The first generalization: Decadic oscillators with q = 2 
and any N 

The decadic anharmonic oscillator exhibits certain solvability features which moti- 
vated its deeper study in non-Hermitian context [18]. The changes of variables make 
this oscillator very closely related to the common quartic problem with a recognized 
relevance of both its non-Hermitian [19] and Hermitian [8, 16, 20] QES constructions. 
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Paying attention to the D ^> 1 domain and abbreviating the parameters s± — s 
and s 2 = t of the respective decadic-oscillator energy and coupling in eq. (30), we 
arrive at the four-diagonal version of our solvability condition (31) at q — 2, 



s 
t 

N — 1 



1 

s 
t 

N-2 



\ 



2 

s 
t 



V 



N- 1 

s 
t 



Po 

Pi 

PN-2 

V Pn-i J 



(38) 



This is the first really nontrivial equation of the class (31). In order to understand 
its algebraic structure in more detail, let us first choose the trivial case with N = 2 
and imagine that the resulting problem 




(with pi ^ due to the definition of N) may be solved by the determination of the 
unknown ratio of the wavefunction coefficients Po/pi = —t from the last line, and by 
the subsequent elimination of t — 1/s using the first line. The insertion of these two 
quantities transforms the remaining middle line into the cubic algebraic equation 
s 3 = 1 with the single real root s — 1. 

The next equation at iV = 3 is still worth mentioning because it shows that the 
strategy accepted in the previous step is not optimal. Indeed, in 



(39) 




(t 2 — s)/2 from the third line 



the same elimination of P\jp2 = —t and of P0/P2 
leads to the apparently ugly result 

st 2 - s 2 - 2t = 0, 
t 3 - 3st + 4 = 0. 

Still, one should not feel discouraged, at least for the following two reasons. Firstly, 
an alternative strategy starting from the elimination of po and P2 leads to the much 
more symmetric pair of the conditions 

t 2 - s 2 t + 2s = 0, 
s 2 -t 2 s + 2t = 
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the respective pre- multiplication of which by t and s gives the difference t 3 = s 3 . 
This means that t — es where the three eligible proportionality constants exist such 
that e 3 = 1. Thus, our problem degenerates to a quadratic equation with the pair of 
the real roots s — t — s^ 1 ' 2 ) such that 

S (D = 2 , S W = -l. 

The second reason for optimism is even stronger: The "ugliness" of the results of 
the elimination may be re-interpreted as an inessential intermediate stage of the 
solution of eq. (39) for all its four unknowns (we may always put Pn = 1) by the 
"brute-force" symbolic manipulations on the computer, by the so called technique 
of Grobner bases [21]. In particular, at N = 3, the computerized experiment of this 
type leads to the step-by-step elimination of the redundant unknowns and to the 
final effective polynomial equation for the single unknown quantity s, 

s 6 - 7s 3 -8 = 0. 

This equation may be verified to possess the same complete set of the real roots as 
above. One may conclude that the real energies of the "strongly spiked" decadic 
oscillator are very easily determined even without a detailed specification of an "op- 
timal" elimination pattern. 

We see that in general one may expect that eq. (38) may give many unphysical 
complex roots as well. This is confirmed by the next step with N = A leading to the 
effective polynomial equation 

s 10 - 27s 7 + 27s 4 - 729s = 0. 

Being tractable by the standard computer software, it results in the set of the mere 
two real roots again, 

s « = 3, s (2) = 0, iV = 4, q = 2. 

It is not difficult to continue along the same path. In general one finds that the 

q = 2 problem may be reduced to a single polynomial equation with 

complex roots s. Still, originally, we were unable to suspect that after all the explicit 
calculations, all the general physical (i.e., real) spectrum of energies proves to be 
quite rich and appears described again by the following closed and still almost trivial 
formula 

s 0) = iV + 2 — 3j, j = 1, 2, ... , j max , j max = entier 

This is our first important conclusion. After one applies the same procedure at the 
higher and higher dimensions N, some more advanced symbolic manipulation tricks 



rn 



i\ + 1 



(40) 
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must be used [21]. Nevertheless, one repeatedly arrives at the confirmation of the 
iV— independent empirical observation (40) and extends it by another rule that at all 
the values of the dimension N, there exist only such real roots that = . This 
means that each " solvability admitting" real energy s requires, purely constructively, 
the choice of its own "solvability admitting" real coupling constant t. 

This is our second important result which parallels completely the similar obser- 
vations made in our preceding paper on the quartic oscillators [20]. Now, a fully 
open question arises in connection with all the q > 2 versions of eq. (31). Do their 
real roots s,t, . . . exhibit the similar pattern as emerged at q = 2? 

4.3 The second generalization: Oscillators with q = 3 and 
their solvability at any N 

At q = 3 we have to solve the five-diagonal eq. (31), 

/ r 1 

s r 2 

t s r 3 

N - 1 t s r 4 

iV-2 t s r 5 

4 t s r 

3 t s 

2 t 

\ 1 

which may be reduced, by means of the similar symbolic computations as above, to 
the single polynomial problem 

t 9 - 12t 5 - 64 i = 
at iV = 3, to the next similar condition 

t 16 - 68 1 12 - 442 1 8 - 501 16 1 4 + 50625 = 
at iV = 4, to the conditions of vanishing of the secular polynomial 
t 25 - 260 1 21 + 7280 1 17 - 1039040 t n - 152089600 1 9 + 2030239744 1 5 + 10485760000 1 
at iV = 5, or to the perceivably longer equation 

t 3 6 _ 777t 3 2 + 135 71 6 t 28 _ 17189 4 6 ot 24 _ 3513570690 1 20 - 

-1198527160446 1 16 + 103857100871252 1 12 + 873415814269404 1 8 + 



N - 1 

r 
s 
t 



( Po 

Pi 

PN-2 

\ Pn-i ) 



(41) 



15 



+74500845455535625 t 4 - 75476916312890625 = 



at N = 6 etc. These computations represent a difficult technical task but at the 
end they reveal again a clear pattern in the structure of the secular polynomials as 
well as in their solutions. One arrives at the similar final closed formulae as above. 
Now one only deals with more variables so that we need two indices to prescribe the 
complete classification scheme 

s = s U) = N + 3 - Aj, 

r = r (i ' fc) = t = t u ' k) = -N - 3 + 2j + 2k, (42) 

k = 1 , 2 , . . . , k max (j ) , k max (j) — N + 2 — 2j , 

- N + 
. 2 

We may re-emphasize that all the real roots share the symmetry r = t but admit 
now a different second root s. The physical meaning of these roots is obvious. Thus, 
the energies of the oscillations in the polynomial well 



j 1? 2, . . . , j max i 3 max CTltiCT 
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will be proportional to the roots r^ ,k \ After the change of variables, the roots 
will represent energies for the alternative, "charged" polynomial potentials 



V^' k=2 \r) = - + ar + br 2 + ... + fr 

r 



etc [15]. 



5 Outlook: QES solutions at q > 4 and selected TV 

5.1 An apparent loss of simplicity at q = 4 and TV < 6 

In our present formulation of the problem (31), we denote the descending diagonals 
as s m with m — 1, 2, 3, 4 and get the equation 
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si 1 




\ 








S2 












S3 


••. — 
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/ Po 


\ 








Pi 






S 4 


Si 










N-1 ■•■ 


s 2 




\ Pn-i 
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S3 
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1 S4 
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Its systematic solution does not parallel completely the above-described procedures. 
In fact, the reduction of the problem to the search for the roots of a single polynomial 
secular equation P(x) = (in the selected auxiliary variable x = — S4) enables us 
only to factorize P(x) on an extension of the domain of integers, 

p(x) = (x + 3) (2 x + 1 - VK) (2 x + 1 + Vh) 

(2 x 2 - 3 x + 3 V5x + 18) (2 x 2 - 3 x - 3 y/bx + 18) 

(2 x 2 - 3 x - v 7 ^ + 8 + 2 v 7 ^) (2 x 2 - 3 x + y/Ex + 8 - 2 Vh) 

(x 2 + x + V5x + 4 + v 7 ^) (x 2 + x- V5x + 4 - v 7 ^) 

(-2^ + 8- 3^ + 3 V5x + 2 x 2 ) (2 ^ + 8- 3^-3 v^rr + 2 x 2 ) 

(-2 v 7 ^ + 8 + 7x - V5x + 2 x 2 ) (2 v 7 ^ + 8 + 7x + v^x + 2 a; 2 ) 

(2x 2 + 2a; + 3- v 7 ^) (2 x 2 + 2 x + 3 + Vt) 

(Vh + 3 - 3 x - y/Ex + 2 x 2 ) (-^ + 3-32; + + 2 x 2 ) 

(2^ + 8-3^ + ^ + 2 x 2 ) (-2^ + 8-33;- v / 5x + 2a; 2 ) . 
From this lengthy formula it follows that we get 

8? = 3, S f = - 1-618, S f = - -0.618 . 

There only exist these three real roots s 4 in this case. 

The similar computerized procedure gave us the real roots also at N = 5. Their 
inspection leads to the conclusion that S2 = S3, si = S4. We did not succeed in an 
application of our algorithms beyond N = 5 yet. Even the N = 5 version of eq. (43) 
in its reduction to the condition 

x 70 - 936 x 65 + 67116 x 60 - 95924361 x 55 - 74979131949 x 50 + 8568894879002 x 45 - 

... - 17459472274501870222336 x 5 + 142630535951654322176 = 

of the vanishing auxiliary polynomial required a fairly long computation for its (still 
closed and compact) symbolic- manipulation factorization summarized in Table 1. 

Marginally, it is worth noticing that the choice of q = 4 is the first instance 
where the popular cubic polynomial forces may emerge via the change of variables 
of ref. [15]. For this reason, in particular, the incomplete character of our q = 4 
solution might prove challenging in the context of the so called PT symmetric quan- 
tum mechanics where the study of cubic oscillators happened to play a particularly 
significant role [22]. 
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5.2 Simplicity re-gained at q = 5 

5.2.1 N = 6 

At q = 5 and N = 6 the symbolic manipulations using the Grobner bases [23] 
generate the secular polynomial in x = s 5 which has the slightly deterring form 

x 91 - 16120 x 85 + 49490694 x 79 - 286066906320 x 73 - 3553475147614293 a; 67 - 

... - 319213100611990814833843025405983064064000000 a; = . 

Fortunately, it proves proportional to the polynomial with the mere equidistant and 
simple real zeros, 

P?\x) =x(x 2 - l) (x 2 - 2 2 ) (x 2 - 3 2 ) (x 2 - A 2 ) (x 2 - 5 2 ) . 

The rest of the secular polynomial is a product of the other two elementary and 
positive definite polynomial factors 

p®( x ) = f[ (x 2 -3kx + 3k 2 ) (x 2 + 3A; 2 ) (x 2 + 3kx + 3k 2 ) 
k=i 

and 

pf = H (x 2 - kx + k 2 ) (x 2 + kx + k 2 ) , 
k=l 

with another positive definite polynomial 

12 

Pi ] = II (x 2 -b k x + c k ) (x 2 + b k x + c k ) 

k=l 

where the structure of the two series of coefficients (see their list in Table 2) is entirely 
enigmatic. 

The subsequent symbolic manipulations reveal a symmetry s 2 = «4 and Si = S5 
of all the real eigenvalues. In the pattern summarized in Table 3, we recognize a 
clear indication of a return to the transparency of the q < 3 results which may be 
written and manipulated in integer arithmetics. 

5.2.2 N = 7 

One should note that in spite of its utterly transparent form, the latter result re- 
quired a fairly long computing time for its derivation. One encounters new technical 
challenges here, which will require a more appropriate treatment in the future [24]. 
Indeed, the comparison of the N = 6 secular polynomial equation with its immediate 
N = 7 descendant 

x 127 - 60 71 x 121 + 1021190617x 115 - 11387407144495a; 109 - . . . + cx ■ 10 6 = 
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shows that the last coefficient 

c = 12537122012272666762007378932665841565459588304127431 1330630729728 

fills now almost the whole line. This case failed to be tractable by our current 
computer code and offers the best illustration of the quick growth of the complexity 
of the q > 5 constructions with the growth of the QES dimension parameter N. 

Fortunately, we are still able to keep the trace of the pattern outlined in Tables 2 
and 3. Indeed, our new secular polynomial factorizes again in the product of four 
factors Pj(x), j = 1, 2, 3, 4 where only the first one has the real zeros, 

p[ 7 \x) = P ( t\x)-(x 2 -Q 2 ) . 
The further three factors fit the structure of their respective predecessors very well, 
pp ( x ) = P^ ( x ) .(x 2 - 9x + 27) (x 2 + 27) (x 2 + 9 x + 27) 

and 

pP = pf . [x 2 - 6 x + 36) (x 2 + 6 x + 36) 

while 

PP = Pf • II (? 2 -f k x + g k )(x 2 + f k x + g k ) . 
k=i 

The subscript-dependence of the new coefficients is listed in Table 4. 

On the basis of the above factorization we may deduce that the pattern of Table 2 
survives, mutatis mutandis, also the transition to iV = 7. Indeed, by inspection of 
Tables 2 and 4 one easily proves that the product function P2(x) Ps(x) P\{x) has no 
real zeros and remains positive on the whole real line of x again. A full parallel with 
the N = 6 pattern is achieved and might be conjectured, on this background, for all 
N, therefore. 

6 Summary 

Our paper offered new closed solutions of Schrodinger equation with polynomial 
potentials at the large angular momenta £ >> 1. This type of construction proves well 
founded and motivated, say, in nuclear physics where, quite naturally, the variational 
calculations in hyperspherical basis lead to the very large values of £ = O(10 3 ) [25]. 
In such a setting, of course, practically any version of the popular \ j£ perturbation 
expansion (a compact review may be found, e.g., in papers [26]) must necessarily 
lead to a satisfactory numerical performance. 

Our present project was more ambitious. We imagined that a rarely mentioned 
week point of all the above perturbative philosophy lies in the notoriously narrow 



19 



menu of the necessary zero-order approximants H [27]. Indeed, in spite of an amaz- 
ing universality of all the different \ji (better known as 1/N) expansion techniques 
(cf. a small sample of the relevant computational tricks in [28]), one usually finds and 
returns to the common harmonic oscillator Hq H °^ = H^ q=0 \ in spite of the wealth 
and variability of the underlying physics [29]. For this reason we recently started 
to study some alternative possibilities offered by the QES models [17, 20]. In our 
present continuation of this effort, a decisive extension of the results of this type is 
given. 

Our text reveals the existence and describes the construction of certain fairly 
large multiplets of "exceptional" £ ^> 1 bound states for a very broad class of poly- 
nomial oscillators. We believe that they might find an immediate application in 
some phenomenological D 1 models where the enhancement of the flexibility of 
the models with q > 1 might lead, say, to a more precise fit of the vibrational spectra 
etc. 

From the mathematical point of view, the most innovative and characteristic 
feature of our new D ^> 1 QES multiplets lies in the existence of the new closed and 
compact formulae for the QES energies and/or couplings at all N. For this reason, 
the corresponding partially solvable polynomial oscillator Hamiltonians H^ q ' N ^ might 
even be understood as lying somewhere in between the QES and ES classes. 

Due to such an exceptional transparency of our constructions of Hq 9 ' N \ a facili- 
tated return to the "more realistic" finite spatial dimensions D = 0(1) might prove 
tractable by perturbation techniques. Two reasons may be given in favor of such a 
strategy. First, due to the specific character of our present "unperturbed" spectra 
and eigenvectors, the perturbation algorithm might be implemented in integer arith- 
metics (i.e., without rounding errors) in a way outlined, preliminarily, in ref. [17] at 
q — 1. Second, the evaluation of the few lowest orders might suffice. This expecta- 
tion follows from the enhanced flexibility of the available zero-order Hamiltonians. A 
priori, a better convergence of the corrections will be achieved via a better guarantee 
of a "sufficient smallness" of the difference between the realistic Hamiltonian H and 
its available approximant H . 

Of course, the detailed practical implementation of the perturbation technique 
represents an independent task which must be deferred to a separate publication. 
With encouraging results, the first steps in this direction have already been performed 
at q = 2 [20]. In parallel, it seems feasible to enlarge further the range of q in zero 
order. Although one has to deal with the fairly complicated symbolic manipulations 
on the computer beyond q = 3, we still intend to perform a deeper analysis of the 
problems with q > 4 in the nearest future [24]. 
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Tables 

Table 1. Sample of the real roots of eq. (43) (q — 4). 
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Table 2. Coefficients b k and c k for q = 5 and N — 6. 
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1 - 3 
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13 


2,5,7 


7-9 
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21 


3,6,9 



Table 3. Real roots of eq. (31) at q — 5 and N = 6 
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Table 4. Additional coefficients f k and at q — 5 and N = 7. 
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1-3 
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28 
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2,8,10 
4,7,11 
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